Molecular Plant

Fruit ripening evolved to be attractive to frugivores that derive energy and nutrition from the fruits in exchange for assisting seed dispersal, which is accompanied by the dramatically change of fruit characteristics, including color, aroma, and texture. The plant hormone ethylene plays a key role in climacteric fruit ripening, while the role of other phytohormones as well as their cross talk with ethylene in modulating fruit ripening remains elusive. Here, we report growth-promoting phytohormone brassinosteroids promote fruit ripening in tomato through regulation of ethylene biosynthesis. Exogenous BR treatment and the increase of endogenous BR content in SlCYP90B3-OE promoted ethylene production and fruit ripening. SlBZR1, a central component and positive regulator of BR signaling pathway, promotes ethylene production and carotenoid accumulation through direct transcriptional regulation of SlACO1, SlACO3 and SlPSY1. Furthermore, SlBIN2, a negative regulator of BR signaling upstream of SlBZR1, decreases ethylene production and carotenoid accumulation. Together, our results demonstrate that BR signaling integrates ethylene and carotenoid biosynthetic pathway to regulate fruit ripening.

Sessile plants exhibit diverse movement behaviors that have long intrigued the scientists. The legume plants display a rhythmic leaflet movement pattern characterized by horizontal opening during the day and vertical closure at night. However, the underlying mechanisms remain largely enigmatic. Here, we isolated a mutant designated as dlm1 (downward leaflet movement1) from Medicago truncatula that displays leaflets downward opening during daytime while upward closure at night. Cellular analyses reveal that this aberrant phenotype correlates with abnormal volume changes of motor cells within the pulvinus. DLM1 exhibits high expression level in the motor organ pulvinus and encodes an unreported nuclear-localized PP2C phosphatase. Notably, phylogenetic analysis demonstrates that species exhibiting rhythmic leaflet movements consistently retain DLM1 homologs, suggesting the possibility of its functionally conserved role in regulating leaflets movement pattern. Structural characterization reveals that DLM1 possesses both phosphatase and kinase domains. Functional complementation assays demonstrate that the phosphatase domain is necessary and sufficient for maintaining the leaflet movement pattern. Collectively, our work uncovers a novel PP2C protein that governs the leaflet movement, providing mechanistic insights into this intriguing plant behavior.

The Asteraceae (Compositae) is the largest flowering plant family, ubiquitous in most terrestrial communities, and morphologically hyper-diverse. An ancient whole genome triplication (paleo-hexaploidization) occurred at approximately the same time as the evolutionary innovation and adaptive radiation of the family during the middle Eocene. Despite its importance, the genomic contents arising from this triplication have yet to be tracked in context of the Asteraceae genome evolution. We applied a synteny oriented phylogenomic analysis of 21 Asterales genomes and to study the paleo-hexaploidization and its consequences to gene, trait, and genome evolution. We identified 15 ancestral linkage groups (ALGs) that date back to the common diploid ancestor of all Asteraceae. Each of these groups was triplicated, resulting in 45 genomic blocks (3x15), which serve as the foundation for cross-family analyses. We demonstrate the complex evolutionary dynamics of the 45 genomic blocks across the Asteraceae phylogeny. We found that modern genomes are genetic mosaics of three progenitor genomes by extensive genomic exchange, chromosomal shuffling and gene fractionation. 157 genes retained three paleo-hexaploid derived syntenic paralogs across most Asteraceae species. Transcription factors (TFs) and auxin-related genes are significantly overrepresented in the conserved triplets, and expression of the paleo-hexaploidy paralogs is spatiotemporally differentiated. These genes are involved in the development of floral capitulum, a remarkable morphological innovation of the family. The discovery of conserved triplicated genes can direct further study to understand the evolutionary innovation, and the synteny-phylogenomic framework and ALGs provide a comparative framework to characterize newly sequenced Asteraceae genomes.

Cross-regulation between hormone signaling pathways is indispensable for plant growth and development. However, the molecular mechanisms by which multiple hormones interact and co-ordinate activity need to be understood. Here, we generated a cross-regulation network explaining how hormone signals are integrated from multiple pathways in etiolated Arabidopsis (Arabidopsis thaliana) seedlings. To do so we comprehensively characterized transcription factor activity during plant hormone responses and reconstructed dynamic transcriptional regulatory models for six hormones; abscisic acid, brassinosteroid, ethylene, jasmonic acid, salicylic acid and strigolactone/karrikin. These models incorporated target data for hundreds of transcription factors and thousands of protein-protein interactions. Each hormone recruited different combinations of transcription factors, a subset of which were shared between hormones. Hub target genes existed within hormone transcriptional networks, exhibiting transcription factor activity themselves. In addition, a group of MITOGEN-ACTIVATED PROTEIN KINASES (MPKs) were identified as potential key points of cross-regulation between multiple hormones. Accordingly, the loss of function of one of these (MPK6) disrupted the global proteome, phosphoproteome and transcriptome during hormone responses. Lastly, we determined that all hormones drive substantial alternative splicing that has distinct effects on the transcriptome compared with differential gene expression, acting in early hormone responses. These results provide a comprehensive understanding of the common features of plant transcriptional regulatory pathways and how cross-regulation between hormones acts upon gene expression.

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in biosynthesis of momilactone A are found in clusters, i.e., MABGCs (Momilactone A biosynthetic gene clusters), in the rice and barnyardgrass genomes. How MABGCs originate and evolve is still not clear. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via horizontal gene transfer (HGT) and followed by recruitment of MAS and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another HGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter being originated from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate HGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.

Plants rely on the innate immune system to sense and respond to a wide range of lifestyle pathogens and to facilitate their survival in natural ecosystems. Pathogen-associated molecular patterns (PAMP)-triggered immunity (PTI) and effector-triggered immunity (ETI) are designated as a two-branched system of innate immunity. Although PTI/ETI share a series of downstream molecular events, systematic analysis of convergent and divergent signaling in PTI/ETI is currently lacking. The phytohormones salicylic acid (SA) and jasmonic acid (JA) are considered to constitute the hormonal backbone of plant immunity, are functionally antagonistic, and play essential roles in defending against biotrophic and necrotrophic pathogens, respectively. However, the distinct performance of two phytohormones in PTI/ETI remains unclear. Here, we systemically investigate and validate the reprogramming of molecular networks during PTI and ETI. Using publicly available Arabidopsis RNA sequence data from 560 samples, we construct a co-expression network under Mock conditions and then explore the differential expression/co-expression changes during PTI, ETI, and Pto DC3000 infection. During PTI and ETI, one-third of genes in the Arabidopsis genome exhibit the same directional differential expression in a manner independent of JA/ethylene/PAD4/SA signaling but show differential co-expression patterns. However, the defense phytohormone network is required for defense against Pto DC3000 infection. We also exhibit the use of this network in prioritizing genes that functioned closely with the proteins directly targeted by elicitors. Overall, this study will deepen our understanding of plant transcriptome in plant immunity and provide new insights into the mode of action of elicitors.

GOLDEN2-LIKE (GLK) transcription factors drive the expression of photosynthesis-associated nuclear genes (PhANGs), indispensable for chloroplast biogenesis. We previously demonstrated that the salicylic acid (SA)-induced SIGMA FACTOR-BINDING PROTEIN 1 (SIB1), a transcription coregulator and positive regulator of cell death, interacts with GLK1 and GLK2 to reinforce their activities. The SIB1-GLK interaction raises the level of light-harvesting antenna proteins in photosystem II, aggravating photoinhibition and singlet oxygen (1O2) burst. 1O2 then contributes to SA-induced cell death via EXECUTER 1 (EX1, 1O2 sensor protein)-mediated retrograde signaling upon reaching a critical level. We now reveal that LESION-SIMULATING DISEASE 1 (LSD1), a transcription coregulator and negative regulator of SA-primed cell death, interacts with GLK1/2 to repress their activities. Consistently, the overexpression of LSD1 represses GLK target genes including PhANGs, whereas the loss of LSD1 increases their expression. Remarkably, LSD1 overexpression inhibits chloroplast biogenesis, resembling the characteristic glk1glk2 double mutant phenotype. The subsequent chromatin immunoprecipitation analysis coupled with quantitative PCR further reveals that LSD1 inhibits the DNA-binding activity of GLK1 towards its target promoters. The SA-induced nuclear-targeted SIB1 appears to counteractively interact with GLK1/2, leading to the activation of EX1-mediated 1O2 signaling. Taken together, we provide a working model that SIB1 and LSD1, mutually exclusive SA-signaling components, antagonistically regulate GLK1/2 to fine-tune the expression of PhANGs, thereby modulating 1O2 homeostasis and related stress responses.

The phytohormone auxin, specifically indole-3-acetic acid (IAA) plays a prominent role in plant development. Cellular auxin concentration is coordinately regulated by auxin synthesis, transport, and inactivation to maintain auxin homeostasis; however, the physiological contribution of auxin inactivation to auxin homeostasis has remained elusive. The GH3 genes encode auxin amino acid conjugating enzymes that perform a central role in auxin inactivation. The chemical inhibition of GH3s in planta is challenging because the inhibition of GH3 enzymes leads to IAA overaccumulation that rapidly induces GH3 expression. Here, we developed a potent GH3 inhibitor, designated as kakeimide (KKI), that selectively targets auxin-conjugating GH3s. Chemical knockdown of the auxin inactivation pathway demonstrates that auxin turnover is very rapid (about 10 min), indicating auxin biosynthesis and inactivation dynamically regulate auxin homeostasis.

Riptortus pedestris (Fabricius) a major soybean pest migrates into soybean fields during pod filling stage resulting in a leaf and stem staygreen while pods without beans syndrome. Given the agricultural importance of this species and the lack of characterized HAMP from piercing-sucking insects we performed a large scale of screening by expression of 87 R. pedestris salivary proteins with signal peptides in Nicotiana benthamiana obtaining a candidate HAMP RPH1. RPH1 activated a series of PTI responses including ROS burst upregulation of defense marker genes such as PR1 WRKY7 WRKY8 Acre31 and CYP71D20 MAPK activation and biosynthesis of phytohormones in plants. RPH1 significantly enhances soybean resistance against R. pedestris feeding. PRR coreceptors BAK1 and SOBIR1 were required for RPH1-induced PTI responses. Remarkably RPH1 homologs were widely distributed in herbivorous insects and majority of homologs from selected species induced cell death or ROS. Thus our results demonstrated that RPH1 is a conserved HAMP within chewing and piercing-sucking insects. We also discovered that R. pedestris evolved four paralogs to overcome the plant immunity triggered by RPH1. This study filled a major gap of HAMP identification from piercing-sucking insect and also deciphered a novel evasion strategy of plant immunity exploited by herbivorous insects. One sentence summaryRiptortus pedestris RPH1, a conserved HAMP in herbivores, activates a variety of PTI responses in plants. To couterdefense, R. pedestris evolved four paralogs to suppress RPH1-induced PTI responses. The author(s) responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (https://academic.oup.com/plcell/pages/General-Instructions) is: Ai Xia (xia@njau.edu.edu).

Receptor-like kinases (RLKs) are essential in plants and phosphorylation is a critical step for their function. Interestingly, RLKs have many non-catalytic kinases/ pseudokinases and the biochemical basis for these pseudokinases remains unclear. FERONIA (FER) is an RLK with kinase activity, but the necessity of its kinase activity for genetic functions has been debated. Here, we uncovered that the kinase-deficient variant FERK565R can activate kinase activity in FER and its homologous through homo/heterodimerization-dependent allosteric activation. We further showed that reactive oxygen species (ROS) significantly promote the dimerization of FER family members. Next, we revealed that mutating the FER P740 within the G-H loop reduces FER dimerization and disrupts its allosteric activation, thus attenuating FERs transphosphorylation for its substrate. This disruption in allosteric activation abolishes the genetic function of FERK565R, impacting ROS production and ABA-mediated stomatal movements. Additionally, we found that MEDOS1 (MDS1), a member of the FER family, is incapable of catalyzing phosphotransfer, but can boost the kinase activity of FER and HERK1 through allosteric activation. These findings settle the debate on FERs inactivated forms, and reveal a new mechanism for allosteric activation of RLKs via redox signaling, enhancing our understanding of pseudokinases in plants. One-sentence summaryFER activates kinase activity of homologous family members through allosteric activation.

Dear Editor, Since the first plant receptor-like kinase (RLK) gene ZmPK1 was cloned from Zea mays in 1990 (Walker & Zhang, 1990), this large gene family has been extensively studied and shown to play crucial roles in growth, development, and immunity (Tang et al., 2017). RLKs are widespread in the plant kingdom, while the biological functions of most RLKs remain largely elusive (Dievart et al., 2020). Given RLKs share a conserved monophyletic RLK/Pelle kinase domain, RLKs in several model plants are classified into distinct families by extracellular domains (ECDs) (Shiu & Bleecker, 2001). However, independent domain shuffling in specific lineages drives the origin of novel families, which raises a question: how about the landscape of RLKs in the whole plant kingdom? Previously, sequence homology-based methods have been widely used for RLK identification and classification, which probably will miss the distantly related proteins but with similar structures and potential novel families unmentioned in the literature. The academic community urgently requires a dedicated database for a systematic overview of the RLK gene family, providing data support for in-depth research on RLK genes. Here, we used a topology-based method to accurately isolate the RLKomes from proteomes. The obtained RLKomes were further classified into (sub)families based on ECD domains. We constructed a comprehensively curated plant RLK database (https://biotec.njau.edu.cn/rlkdb/), which contains valuable resources for investigating the origin and evolution of the RLK family and multiple online tools for personalized analysis.

The red/far-red sensing photoreceptor phytochrome B (phyB) governs multifaceted plant development and responses to light and temperature stimuli. PhyB photoconversion between red-absorbing, inactive Pr and far red-absorbing, active Pfr states, imparted by its covalently bound bilin chromophore, enables rapid switching and plasticity of phyB signaling activities. The phyBY276H variant (YHB) is photochemically inert but adopts a constitutively active Pfr-like structure regardless of light conditions, which becomes a versatile model to dissect phyB signaling mechanisms. Here, we conducted a large-scale EMS mutagenesis screen on YHB-expressing transgenic lines, mining intragenic suppressor mutations that would unveil critical residues for phyB structure-function relationships. Comparative analyses of 26 nonsense variants suggested modular organization of phyB overall structure and dispensability of the C-terminal HKRD domain for phyB signaling. Amongst fourteen novel and nine known loss-of-function missense variants identified herein, G284E was of particular interest for its fully suppressed constitutive activity in darkness and its restored photochemistry and light responsiveness. The G284E mutation was further tested to also nullify another constitutively active phyBY303V allele by eliminating chromophore attachment. P309L was the sole variant identified which fully suppressed YHB in both dark and light conditions. C402Y profoundly elicited YHB protein instability. Three variants G118R, C402Y and G538D markedly reduced chromophorylation levels of YHB. Although the chromophore binding site variant C357Y was a strong loss-of-function allele, it retained residual signaling activity with respect to PIF3 protein turnover in dark-grown seedlings, presumably due to its ability to noncovalently bind chromophore. Two tandem prolines (P799, P800) proved critical to YHB structural integrity/stability as well as signaling activity. In summary, these diverse variants shed new insights into multiple levels by which the YHB (and thereby phyB) signaling is initiated, tuned, and disseminated.

Auxin and auxin-mediated signaling pathways involved in the regulation of lateral root development are well documented. Although exocytic vesicle trafficking plays an important role in auxin efflux carriers PIN recycling, and polar auxin transport during lateral root formation, however, the mechanistic details of these processes are not well understood. Here, we demonstrate that BYPASS1-LIKE (B1L) regulate lateral root development via exocytic vesicular trafficking-mediated polar auxin transport in Arabidopsis. In b1l mutants, the number of lateral roots increased significantly, and the phenotypes were mainly attributed to lateral root primordium initiation but not to the defects in lateral root primordium development. Furthermore, the auxin signal was stronger in the lateral root primordium of the b1l mutant at stage I than those observed in the wild-type (WT). Moreover, exogenous auxin and auxin transport inhibitory treatments indicated that the phenotype of lateral roots in b1l mutants can be attributed to higher auxin levels and that B1L regulates auxin efflux. Consistently, auxin efflux carriers PIN1-GFP and PIN3-GFP were expressed at higher levels in the lateral root primordium of the b1l mutants. Interestingly, we found that B1L interacted with the exocyst and b1l mutant showed a defect in PIN2 exocytosis. Finally, we found that B1L cooperated with EXO70B1 to regulate lateral root formation. Our findings reveal an essential regulatory mechanism of B1L that interacts with the exocyst to regulate PIN-mediated polar auxin transport and lateral root initiation.

Pyridoxal-5-phosphate (PLP), the enzymatic cofactor form of Vitamin B6 (vitB6), is a versatile compound that has essential roles in metabolism. Cellular PLP homeostasic regulation is currently not well understood. Here we report that in Arabidopsis, biosynthesized PLP is sequestered by specific aminotransferases (ATs), and that the proteins ROOT UV-B SENSITIVE 1 (RUS1) and RUS2 function with ATs to regulate PLP homeostasis. The stunted growth phenotypes of rus1 and rus2 mutants were previously shown to be rescuable by exogenously supplied vitB6. Specific residue changes near the PLP-binding pocket in ASPARTATE AMINOTRANSFERASE2 (ASP2) also rescued rus1 and rus2 phenotypes. In this study, saturated suppressor screens identified 14 additional suppressor of rus (sor) alleles in four aminotransferase genes (ASP1, ASP2, ASP3, or ALANIN AMINOTRANSFERASE1 (AAT1)), which suppressed the rus phenotypes to varying degrees. Each of the sor mutations altered an amino acid in the PLP-binding pocket of the protein, and sor proteins were found to have reduced levels of PLP conjugation. Genetic data revealed that the availability of PLP normally requires both RUS1 and RUS2, and that increasing the number of sor mutants additively enhanced the suppression of rus phenotypes. Biochemical results showed that RUS1 and RUS2 physically interacted with ATs. Our studies suggest a mechanism in which RUS1, RUS2 and specific ATs work together to regulate PLP homeostasis in Arabidopsis.

Salicylic acid (SA), a central hormone in plant immunity, is biosynthesized via a recently elucidated phenylalanine-derived pathway in most seed plants. This pathway requires benzyl alcohol as a key substrate for the formation of the SA precursor benzyl benzoate. However, how benzyl alcohol is produced in plants was unclear. Here, we identify a two-step conversion of benzoyl-CoA to benzyl alcohol via benzaldehyde in Nicotiana (N.) benthamiana. From a forward genetic screen for SA-deficient mutants, the and {beta} subunits of heterodimeric benzaldehyde synthase (BalS) involved in the conversion of benzoyl-CoA to benzaldehyde were found to be required for SA biosynthesis in N. benthamiana. Further reverse genetic analysis revealed that the NADPH-dependent benzaldehyde reductase (BalR1) acts downstream of BalS to convert benzaldehyde to benzyl alcohol. Interestingly, OsBalR1, but not OsBalS or OsBalS{beta}, is required for maintaining high basal SA levels in rice, suggesting the presence of redundant benzoyl-CoA-reducing activities or alternative biosynthesis routes for benzyl alcohol production. Together, this work defines the missing enzymatic steps in phenylalanine-derived SA pathway and provides insights into the evolutionary diversification of SA production strategies in plants.

Brassinosteroids (BR) and Target of Rapamycin Complex (TORC) are two major processes coordinating plant growth and stress responses. BRs function through a signaling pathway to extensively regulate gene expression and TORC is known to regulate translation and autophagy. Recent studies revealed that these two pathways crosstalk, but a system-wide view of their interplay is still missing. Thus, we performed transcriptome, proteome, and phosphoproteome profiling of Arabidopsis mutants with altered levels of either BIN2 or RAPTOR1B, two key players in BR and TORC signaling, respectively. We found that perturbation of BIN2 or RAPTOR1B levels affects a common set of gene-products involved in growth and stress responses. Additionally, we performed Multiplexed Assay for Kinase Specificity (MAKS), which provided a system-wide view of direct BIN2 substrates. Furthermore, phosphoproteomic data was used to reconstruct a kinase-signaling network and to identify novel proteins dependent on BR and/or TORC signaling pathways. Loss of function mutants of many of these proteins led to an altered BR response and/or modulated autophagy activity. Altogether, these results provide genome-wide evidence for crosstalk between BR and TORC signaling and established a kinase signaling network that defines the molecular mechanisms of BR and TORC interactions in the regulation of plant growth/stress balance.

The levels of the important plant growth regulator indole-3-acetic acid (IAA) are tightly controlled within plant tissues to spatiotemporally orchestrate concentration gradients that drive plant growth and development. Metabolic inactivation of bioactive IAA is known to participate in the modulation of IAA maxima and minima. IAA can be irreversibly inactivated by oxidation and conjugation to Aspartate and Glutamate. Usually overlooked because its reversible nature, the most abundant inactive IAA form is the IAA-glucose (IAA-glc) conjugate. Glycosylation of IAA is reported to be carried out by the UDP-glycosyltransferase 84B1 (UGT84B1), while UGT74D1 has been implicated in the glycosylation of the irreversibly formed IAA catabolite oxIAA. Here we demonstrate that both UGT84B1 and UGT74D1 modulate IAA levels throughout plant development by dual IAA and oxIAA glycosylation. Moreover, we identify a novel UGT subfamily whose members modulate IAA homeostasis during skotomorphogenesis by redundantly mediating the glycosylation of oxIAA.

Stomata play critical roles in gas exchange and immunity to pathogens. While many genes regulating early stomatal development up to the production of young guard cells (GCs) have been described in Arabidopsis, much less is known about how young GCs develop into mature functional stomata. Here we performed a maturomics study on stomata, with "maturomics" defined as omics analysis of the maturation process of a tissue or organ. We developed an integrative scheme to analyze three public stomata-related single-cell RNA-seq datasets and identified a list of 586 genes that were specifically up-regulated in all three datasets during stomata maturation and function formation. The list, termed sc_586, is enriched with known regulators of stomatal maturation and functions. We selected two candidate G2-like TFs genes, MYS1 and MYS2, from the list to investigate their roles in stomata. Our results showed that these two genes redundantly regulate the size and hoop rigidity of mature GCs, and their double mutations caused mature GCs to have severe defects in regulating their stomatal apertures. Our analysis thus provides a valuable gene list for studying GC maturation and function formation.

Carotenoids are health-promoting plastidial isoprenoids with essential functions in plants as photoprotectants and photosynthetic pigments in chloroplasts. They also accumulate in specialized plastids named chromoplasts, providing color to non-photosynthetic tissues such as flower petals and ripe fruit. Carotenoid accumulation in chromoplast requires specialized structures and proteins such as fibrillins. Although fibrillins were first reported as structural components of carotenoid sequestering structures in chromoplasts, later work revealed roles in chloroplasts and other plastid types. However, the association of fibrillins with carotenoids in plastids other than chromoplasts has remained unexplored. Here we show that a member of the fibrillin family, FBN6, interacts with phytoene synthase (PSY, the first committed and rate-determining step of the carotenoid pathway) to promote its enzymatic activity. Transient overexpression of FBN6 in Nicotiana benthamiana leaves results in a higher production of phytoene, the product of PSY activity, whereas loss of FBN6 activity in Arabidopsis thaliana mutants dramatically reduces the production of carotenoids during seedling deetiolation and after exposure to high light. Our work hence demonstrates that fibrillins not only promote the accumulation of carotenoids but also their biosynthesis.

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species rarely occur in nature mainly due to reproductive barriers but how such hurdles can be overcome is largely unknown. xBrassicoraphanus is a fertile intergeneric allopolyploid synthesized between Brassica rapa and Raphanus sativus in the Brassicaceae family. Genomes of B. rapa and R. sativus are diverged enough to suppress synapsis formation between non-homologous progenitor chromosomes during meiosis, and we found that both genomes reside in the single nucleus of xBrassicoraphanus without genome loss or rearrangement. Expressions of syntenic orthologs identified in B. rapa and R. sativus were adjusted to a hybrid nuclear environment of xBrassicoraphanus, which necessitates reconfiguration of transcription network by rewiring cis-trans interactions. B. rapa coding sequences have a higher level of gene-body methylation than R. sativus, and such methylation asymmetry is maintained in xBrassicoraphanus. B. rapa-originated transposable elements were transcriptionally silenced in xBrassicoraphanus, rendered by gain of CHG methylation in trans via small RNAs derived from the same sequences of R. sativus subgenome. Our work proposes that not only transcription compatibility but also a certain extent of genome divergence supports hybrid genome stabilization, which may explain great diversification and expansion of angiosperms during evolution.